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ABSTRACT 

A signihcant percentage of OB stars are runaways, so we should expect a 
similar percentage of their evolved descendants to also be runaways. However, 
recognizing such stars presents its own set of challenges, as these older, more 
evolved stars will have drifted further from their birthplace, and thus their veloc¬ 
ities might not be obviously peculiar. Several Galactic red supergiants (RSGs) 
have been described as likely runaways, based upon the existence of bow shocks, 
including Betelgeuse. Here we announce the discovery of a runaway RSG in M31, 
based upon a 300 km s“^ discrepancy with M31’s kinematics. The star is found 
about 21'(4.6 kpc) from the plane of the disk, but this separation is consistent 
with its velocity and likely age (~10 Myr). The star, J004330.06-1-405258.4, is an 
M2 I, with My = —5.7, logL/L0=4.76, an effective temperature of 3700 K, and 
an inferred mass of 12-15Mq. The star may be a high-mass analog of the hyper¬ 
velocity stars, given that its peculiar space velocity is probably 400-450 km s“^, 
comparable to the escape speed from M31’s disk. 
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Introduction 


Massive stars form in giant molecular clouds, creating OB associations (e.g. JLada fc Lada 


2003|), the members of which share similar space velocities. OB stars that are spatia lly close 
to one of these associations, but with discrepant radial velocities, were hrst noted bv iBlaauw 
(1l96lh . who termed these “runaways” and proposed a possible explanation for their origin: 
if the primary in a binary system sheds a significant fraction of its mass (say, in a SN ex¬ 
plosion), then the secondary would be set loose with nearly its orbital velocity. This would 
cause the star’s radial velocity to disagree with that of other members of the association, 
and over time, would cause the star to move awa y from its fellow’ s . A ra dial velocity study of 
a large sample of Galactic runaway OB stars by iGies fc Bolton! fjl986h effectively ruled out 
this explanation, and instead favored a dynam ical evolution explanation, a sce nari o that has 


been r ecently supported by the simulations of iFnjii fc Portegies ZwartI (120111) and lOh et ah 
(l2ni5[) . 


Regardless of thei r origins, a sigiiificaii t percentage of OB stars are considered runaways 


(10-50% according to iGies fc Boltonlll986l. and references therein). Yet, few evolved mas¬ 


sive stars have ever been identified as runaways, presumably as these stars are older, and 
as a star moves farther away, it begins to lose the context of its birth association. With¬ 
out that context, unusual velocities become harder to distinguish. Nevertheless, a few red 


(Noriesra-GresDO et ah 

1997a: 

Mackev et al 

. 2012) 

some OB runaways ( 

Norieea-GresDO et ah ! 

997hl). 


a star is a runa way. Galactic RSGs con sidered to be runaways include Betelgeuse, fi Gep, 
and IRG-10414 f Gvaramadze et ah 2014 ). 


While conducting a radial velocity study of RSGs in M31, we discovered that the RSG 
J004330.06-1-405258.4 is not only quite isolated, but also possesses a radial velocity 300 km s“^ 
at odds with that expected from its location. We conclude that this is a runaway RSG, the 
first identified in another galaxy, and the fastest known. 


2. Observations 


J004330.06-1-4 05258.4 was one of many M31 RSGs we observed spectroscopically. Our 
sample came from iMassey et ahl (120091) . who had identified candidate RSGs a nd foreground 


stars in M31 using the photometry of the Local Group Galaxy Survey iLGGS. lMassey et ah 


20061 ). For very red stars, V-R remains a temperature discernment, but fi-Fbecomes pri- 
marily an indicator of surface gravity due to the effects of line-blanketing in the B bandpass 
(iMassey et al.ll2006l) . Thus, red supergiants can be separated from foreground stars in a 
B-V, V-R diagram. Gandidate RSGs were separated from the candidate foreground stars 
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using the following equation: 

B-V = -1.599(1/ - Rf + 4.18(1/ - R) - 0.83 (1) 

We summarized the available photometry for J004330.06+405258.4 in Tabled! The location 
of J004330.06+405258.4 in the two-color diagram of our sample is shown in Figured! where 
we have colored the expected RSGs in red and the expected foreground stars in black, based 
upon the above cut-off. The location of our runaway star, J004330.06+405258.4 is very much 
in the expected RSG area. 


The full results of the radial velocity study will be discussed elsewhere, but here we 
will summa rize the data relevant to J004330.06+405258.4. The object was observed using 
Hectospec fiFabricant et al.l 120051) . a 300 optical hber-fed spectrograph on the 6.5-m MMT 
telescope. We used the 270 line mm“^ grating, which is blazed at 5000 A, and covers a 
spectral range from 3650 to 9200 A. The grating provides a dispersion of 1.2 A pixeR^ and 
a spectral resolution of 6.2 A. Some of the hbers were assigned to blank sky to be used for 
sky subtraction, and calibration included HeNeAr and quartz lamp exposures. Our runaway 
was observed during four different nights in the Fall of 2014 in queue mode for 90 minutes 
of exposure each, as summarized in Table [2l Following the observations, data were passed 
through the SAO pipeline, and the wavelength zero-points adjusted slightly by using night 
sky lines. The wavelengths were then corrected to a heliocentric reference frame. 


As no blocking hlter could be used, we expect contamination by second-order blue light 
at the longer wavelengths of our spectrum. However, at the Ga ll triplet lines (8498, 8542, 
8662 A), which we use for radial velocities, we expect this contamination will be only ~3%, 
as B — R 2.9 (Tabled]), implying that the overlapping second-order blue at 4250 A will 
be roughly a factor of 15 smaller in flux at 8500 A, and the light dispersed by an additional 
factor of 2 in second order. 


Flux calibration was provided by observations of the spectrophotometric standard Feige 
34 midway through the semester. The resulting sensitivity curves were kindly provided by 
Nelson Galdwell. Our experience is that this is typically good to 5-10% in the fluxes relative 
to wavel engths (i.e., the color s), although occasionally much larger errors may be present 
(see also IFabricant et al.l 120081) . 


3. Analysis 

The hrst s tep in our larger pro ject is to establish whether the stars selected from the 
photometry by iMassev et al.l (120091) were truly RSGs in M31 or merely foreground stars. 
Except for stars in the NE section of the galaxy, the combination of the '^^250 km s“^ 
rotational velocity and 300 km s“^ systemic velocity (jde Vaucouleurs et al.lll99ll) makes 
it easy to separate fore ground Milky Way dwarfs from bona fide M31 members using radial 
velocities, as shown by iDrout et ahl (120091) for yellow supergiants. 
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The observed radial velocities were calculated using XCSAO, a cross correlation tool, 
in irafQ. For our velocity templates, we used 21 Hectospec s pectra of six M31 RSGs for 
which radial velocities were already known flMassey et alJl2nn9h . and we restricted the cross 
correlation to the Ca ll triplet lines, which are very strong in RSGs; this avoids using the 
extremely broad molecular features. We then averaged the result for each spectra, weighting 
the results in accord with the the internal errors. For stars with more than one observation, 
the results of each spectra were averaged, with weights assigned by the weighted errors of 
each spectra. Full details will be given in the larger paper that is in preparation. 


We then compared the observed radial velocities with tha t expected f r om th e star’s 
location in M31 using the simple ki nematic model adopt ed by iDront et al.l (120091 ). which 
is based upon the seminal study bv iRubin fc FordI fjl970l) . In general, the expected radial 
velocity Wxp from circular rotation will be Kxp = Vq + R(i?) sin,^cos6', where Vq is the 
systemic radial velocity, V(R) is the rotational velocity at a distance from the galactic center 
(R) within the plane of the disk, ^ is the angle between the line of sight and the perpendicular 
to the plane of the galaxy, and 9 is the angle from the semi-major axis. Thus, cos6' is X/R, 
where X is the position along the major axis. The gross simplihcation that V(R) is a constant 
('whi ch is equivalent to s aying that dark ma tter dominates the kinematics) works remarkably 
well (iDrout et al.l 120091: iMassev et al.l 1200911 . leaving us then wi th a linear relat i onshi p. For 
consistency with these earlier works, and to simplify using the iRubin fc FordI (1197011 data, 
we have adopted ^=77° and then used Wxp = —295 -|- 241.5(X/i?). This eq uation does an 
excelle nt job of htting both the velocities of the Hll re gions measured by 


Rubin & Ford 


()l970l) as well as the RSG radial velocities measured by iMassev et al.l (120091) . as shown in 
Figure 2 of the latter work. 


We compare the velocity difference (Wbs — F)jxp) for each star in Figure O The suspected 
RSGs and the suspected foreground stars split generally well into two distinct groups. The 
RSGs cluster around a velocity difference of zero, indicating that their observed velocities 
are about equal to the velocities that we expect based on M31’s kinematics. The foreground 
stars have velocity differences that are about equal to the negative of Wxp, as their Wbs 
values are about zero. It is immediately clear that one star stands out as having a very 
peculiar negative velocity. This star is J004330.06-1-405258.4, which we have highlighted 
with a green pentagon. There are a few other outliers: three purported RSGs that are 
obviously foreground stars (due to being marginal in the 2-color diagram), and several stars 
whose photometry was compromised by crowding and/or whose spectra were poor. These 
will be discussed in more detail in our larger paper. 


The radial velocities of J004330.06-1-405258.4 are given in Tabled and we can see that 


^IRAF is distributed by the National Optical Astronomy Observatory, which is operated by the Associ¬ 
ation of Universities for Research in Astronomy (AURA) under a cooperative agreement with the National 
Science Foundation. 
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the four measurements are quite consistent. Furthermore, the spectra of this star are well 
exposed, and the triplet lines readily measured; see Figure 01 The observed radial velocity 
of J004330.06+405258.4 is —630 km s“^, while its expected radial velocity is —328 km s“^, 
leading to a —302 km s“^ peculiar radial velocity. 


Consistent with its peculiar radial velocity, J004330.06+405258.4 is also very well sepa¬ 
rated spatially from other massive stars, as shown in Figure [3] from the LGGS. We measure 
a separation of about 21' (4.6 kpc) from the major axis. Is this separation reasonable? 
Assuming an age of 10 Myr, and a tangental velocity similar to the peculiar component of 
its radial velocity (i.e., —300 km s“^), we can expect J004330.06+405258.4 to have traveled 
~3 kpc from its birthplace, an extremely good match to what we observ^. 


3.1. Physical Properties of the Runaway 


We assigned a spectral type of M2 I to J004330.06+405258.4. M-type subclasses are 
determined primarily by the depth of the TiO bands. Rather than normalize the spectra 
(a problematic issue for very red stars where there is only “pseudo” continua), we instead 
coi npared the spec t ra of J004330.06+4052858.4 to stars that had been previously classified 
bv iLevesque et ahl (120051) digitally in log space. In Figure [5] we show that the spectrum is 
well matched to that of M2 I star BD+59°38. 

We can determine the effective temperat ure Tgfj of J004330.06+405258.4 by using the 
MARG S stellar models (IGustafsson et al.ll2008l ). following the same procedure as in iMassev et al 
(120091) . We adopted a log^f of 0.0, and compared the spectral features (primarily the depths 
of the molecular bands) to those of models of different effective temperatures, adjusting the 
color excess as needed. A “good” fit was determined using the TiO bands between about 
6000 and 7000A. An example is shown in Figure [6l The results of these fits are shown in 
Ta ble El The ternperat ure of 3700 K can be compared to the 3675 K temperature found 
by iMassev et al.l (120091 1 as the median of M2 I stars in M31 using the same super-solar 
metallicity models. 

With this value for the effective temperature we can now determine the bolometric 
luminosity. We measured E{B — V) values of 0.15-0.20 from our htting. This amount of 


reddening is consistent with the typical 0.13 value found from OB stars by iMassev et al. 


(120071) . but given the star’s location far from any star-forming region, i s probably indicativ e 


of some circumstellar component, as is commonly found with RSGs (IMassev et al 


2005 1. 


Adopting Ay = 3.1 x 0.15 = 0.46, and a distance of 760 kpc ffrom Ivan den BerghI 1200011 . 
we hnd an absolute visual magnitude My = —5.65. The bolometric correction at V for 


^It takes 11 Myr for a IS Mf?^ star to beiiiK co re He-burning, and 15 Myr for a 12Mq according to the 
Geneva evolutionary models (jEkstrdm et alJl2012r) 
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the 3700 K MARCS model is —1.51, giving us a bolometric magnitude Mboi = —7.16, or 
logL/LQ of 4.76. Alternatively, we can use the 2MASS phot ometry (given in Ta ble [H) to 
hnd the bolometric luminosity. We expect Ak ~ 0.112 x Ay fjSchlegel et ahlllOOSl ) or 0.05. 
Thus, Mk = —9.93. The bolometric correction at K for a star of 3700 K is +2.76 according 
to the MARCS models. Thus from the K-band, we derive Mboi = —7.17, in near perfect 
agreement with the value derived from V. A comparison with the evolutionary tracks of 


Ekstrom et al.l (12012l ) gives an inferred mass of J004330.06+405258.4 of 12-15M0. 


4. Discussion 


We have established that J004330.06+405258.4 is runaway RSG, the hrst such star 
clearly identihed in another galaxy. Furthermore, with a peculiar radial velocity of-300 km s“^ 
it is the fastest known runaway massive star of which we are aware. Given its distance, a 
direct measurement of its tangental component via proper motions is not practical given 
current methods, but the spatial separation from M31’s disk (Figure [3]) suggests that the 
tangential component of the velocity is similar, and thus that the peculiar space velocity is 
400-450 km s'h 


Lower mass “hypervelocity” stars are known in our own Galaxy; these are stars that 
are moving at peculiar velocities of >500 km s“^, and, beca use of their rnuch g reater ages. 


have traveled many tens of kpc in their lifetimes (see, e.g., iBrown et al.ll2007(l . They are 


escaping from the Galaxy. Our M31 RSG runaway may be a high-mass analog of such stars, 
rather than related to traditional runaways, given that its vel ocity is so mu ch greater than 
the 30 km s“^ that usually use to distinguish an OB runaway (jBlaauwlll96ll ). We note that 
the peculiar velocities o f the other known Galac tic RSGs are quite m odest by compa rison: 
IRG-10414 at 70 km s~^ JGvaramadze et ah 2014), u C ep at 22 km s~^f Cox et ah 2012), and 
Betelgeuse at 56 km s~^(Noriega-GresDo et al.Tl997a h Similarly only two 0-type stars in 


the classic ICruz-Gonzalez et a1 
and both are under 120 km s“^ 


f 


(119741) catalog have peculiar radial velocities >100 km s 


-1 


We were curious what evidence there was of si milar objects i n other galaxies. In their 
study of yellow supergiants (YSGs) in the SMC, iNengent et all (120101 ) do, in fact, note 
that one star, J01020100-7122208, has an anomalously high radial velocity, +307 km s“^, 
compared to the SMCs systemic velocity of 158 km s“^. They suggest the star is either a 
binary or a runaway; based on a single observation, they cannot tell, but it is clear from 
their Figure 7 that the stars radial velocity stands out by 100 km s“^ from the other YSGs. 
Similarly, the histograms of the radial velocities of yellow and red supergiants in the LMC 
studied bv iNengent et ahl (120121 ) also show two interesting outliers (their Figure 4): there is a 
YSG, J04530398-6937285, that has a radial velocity of +373 km s“^, and a RSG, J04482407- 
7104012 (CPD-71°285), with a radial velocity of +401 km s“^, both of which can be compared 
to the LMCs systemic velocity of +278 km s“^. Again, further observations would be needed 
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to establish that these velocities are not high due to binary motion. Amongst the RSGs in 
M33, there is an intriguing example, J01 3403.34+302611.7, with a radial velocity that is 


discrepant by 160 km s ^ (see Figure 8 in iDront et al.ll2ni2h . None of these are, of course. 


as extreme as the case of our M31 RSG runaway. 

Gould J004330.06+405258.4 be a lower mass object? Any alternative explanation needs 
to account not only for the stars radial velocity but also its photometry. We can essentially 
rule out the star being a foreground dwarf: an M2 Galactic dwarf at V = 19.2 would have 
to be at 0.7 kpc to be this faint, and with a radial velocity of —630 km s“^ it would have 
to be a hypervelocity star; an unlikely coincidence given how close it is in the plane of the 
sky to M31’s disk. (We would also not expect it to be as reddened.) What if it were an 
asymptotic giant branch (AGB) star in M31’s halo? The bolometric magnitude Mboi = —7.2 
(logL/L 0 = 4.76) is faint enough that we expect some contamination by intermediate-mass 
AGBs in this magnitude range. However , the typical veloc ity dispersion of halo stars in 
our own Milky Way is about 100 km s“^ flBrown et al.ll2010l) . and we expect that M31’s is 
not much different. We would then still be left trying to explain why the star had a radial 
velocity 330 km s“^ more negative than M31’s —300 km s“^ systemic velocity, as this is far 
in excess than what we expect for a typical halo object. 

J004330.06+405258.4 is seen in relative close proximity (in projection) to the galaxy 
M32, the compact low-luminosity elliptical visible ~ 8 ' to the west in Figure [21 Gould there 
be a connection? It is well known that dwarf ellipticals have very little or no current star 
formation. In th e case of M32, about 3% of its population consists of “young” (ages <2 Gyr), 
metal-rich stars flMonachesi et al.ll2012l) . but there is no evidence of any massive stars. Nor 
is it likely that J004330. 06+405258.4 is an AGB ass ociated with M32, as the galaxy’s radial 
velocity is -205 km s“^ (jde Vauconlenrs et al.lll99lh . which would leave us an even greater 
discrepancy in radial velocity. 

Is J004330.06+405258.4 escaping from M31? This question is largely moot, as the star 
will only live another million years or so before undergoing a SN explosion; in that time it 
will only move only another 400-450 pc. We do note, however, that it is moving sufficiently 
fast to have escaped the gravitational attraction of the disk. Although the gravitational well 
of M31 is complex, with a i nassive dark-rnatter h alo, we can make a crude estimate simply 
using the mass of the disk. [Rubin &: FordI fll970[) calculates that the mass within 24 kpc of 
the center of M31 is ~ 2 x 10^^ Mq. The projected distance of J004330.06+405258.4 is 5.6 
kpc. We would thus expect naively the escape velocity is on the order oi 560 km s“^. (We 
note as a rea lity check that the escape velocity for the sun from the Milky Way is 550 km s“^ 
according to Smith et al. 2007.) Thus it is not unreasonable that J004330.06+405258.4 has 
gotten to where it is today. 


The presence of bow shocks in association with runaway Galactic RSGs (iMackev et al. 


20121 ) leaves a question as t o whether or not we could de tect such a structure around our 
M31 runaway. According to [Noriega-Grespo et al.l fjl997al) . Betelgeuse has a bow shock 0.8 

























pc in size. At the distance of M31, a similiar-sized bow shock would extend about 0.2". This 
size is large enough to be readily detectable, at least using space-based imaging. However, it 
remains to be seen what the effect of our runaway’s increased speed and the likely decreased 
hydrogen column density would be on bow shock size and presence. 
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Table 1. Photometry of J004330.06+405258.4 


Measurement Value 

(mags) 

Error 

(mags) 

pa 

19.212 

0.005 

B-V^ 

1.933 

0.013 

V-R^ 

1.077 

0.006 

iV 

14.525 

0.087 

J-I& 

1.076 

0.106 


^Fr om 

2006 h . 


the LGGS flMassev et ah 


’’From 


2MASS f Gutri et al.ll2003l) 


Table 2. Radial Velocities and Results of the Model Fitting 


HJD 

Radial Velocity (km s ^) 

MARGS htting 

Kbs 

a 


E{B - V) 

2456926.254 

-625.9 

1.02 

3725 K 

0.15 

2456982.310 

-633.7 

1.12 

3700 K 

0.20 

2456987.260 

-623.6 

1.73 

3700 K 

0.15 

2456989.177 

-632.7 

1.01 

3700 K 

0.20 

Adopted Values 

-630 

1 

3700 K 

0.15 
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0.8 1 1.2 1.4 1.6 1.8 

V-R 


Fig. 1.— Two-color diagram of red stars observed with Hectospec. The objects in red 
are suspected red supergiants, and the objects in black are suspected foreground stars. 
The assumed division between the two is given by Equation [H J004330.06-1-405258.6 is 
represented by the green pentagon. 
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Fig. 2.— Difference between the observed velocity (14bs) and the expected velocity (14xp)- 
The stars expected to be RSGs on the basis of their location in the two-color diagram 
(Fignre[T]) are shown with red symbols; as expected, these cluster around a velocity difference 
of 0. The stars expected to be foreground stars based upon their photometry are shown as 
black symbols, and as as expected these cluster around around a line with slope -1 in this 
diagram. The runaway star is noted with a green pentagon. 
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Fig. 3.— Location of J004330.06+405258.4 in M31. Here J004330.06+405258.4 is high¬ 
lighted with a blue circle. The distance from the runaway to the semi-major axis is about 
21' (4.6 kpc). This image was taken in the H-band as part of the LGGS. N is up and E is 
to the left. 
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Fig. 4.— Spectrum of J004330.06+405258.4. The Ca ll triplet used for the radial velocities 
are clearly visible at 8498-8662A. 
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Fig. 5.— Comparison with an M2 I Spectral Standard. A spectrum of J004330.06+405258.4 
('bl ack') is compared with a spectrum of BD+59° 38 (red). The latter was classihed as M2 I 
bv lhevesque et al.l (120051) . 
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Fig. 6.— Model Fit. A spectru m of J004330.06+405 258.4 (black) is compared with a 3700 K, 
log g=0.0 MARCS model from Massev et al. (2009). 












